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Abstract 

We present the simulations of heavy-ion collisions in terms of participant- 
spectator matter. We find that this matter depends crucially on the 
collision dynamics and history of the nucleons. The important changes 
in the momentum space are due to the binary nucleon-nucleon colli- 
sions experienced during the high dense phase. This was otherwise 
not possible with mean field alone. The collisions push the colliding 
nucleons into midrapidity region responsible for the formation of par- 
ticipant matter. This ultimately leads to thermalization in heavy-ion 
collisions. 
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1 Introduction 



The heavy-ion collisions at intermediate energies are the center of present day 
nuclear research. This is because of several rare phenomena emerging at these 
energies and also their utility in several other branches of physics. [TJ, El [3] The 
study of hot and dense matter and its relation to nuclear equation of state 
and cross section has always fascinated the researchers. [TJ El [3J HI |6] Some 
of the rare phenomena emerging are the multifragmentation, [TJ El [3] collec- 
tive flow, [21 [3] stopping^ [3] as well as subthreshold particle production!?] 
etc. In the low energy heavy-ion collisions, the fusion and decay of excited 
compound nucleus and fission dominate the physics. jS] Whereas at high inci- 
dent energies, the transparency and complete disassembly of nuclear matter 
happens. In terms of theoretical description, the relative dominance of real 
and imaginary parts of G- matrix decides the fate of a reaction. [Qj At low 
incident energies, the Pauli principle forbids the nucleon-nucleon binary col- 
lisions and scattering. [TJ El E] The attractive mean field decides the fate of 
colliding pairs. This results into fusion and consequently, into its decay. The 
frequent nucleon-nucleon collisions at high incident energies make imaginary 
part very significant. However, both real and imaginary parts of complex G- 
matrix are equally important at intermediate energies. This picture can also 
be looked in terms of participant-spectator matter and fireball concept. [U [TO] 
The nucleon-nucleon collisions play an important role in destroying the mu- 
tual initial nucleon-nucleon correlations and memory of nucleons. This ulti- 
mately leads to the above rare phenomena at intermediate energies. 

We here plan to understand the formation of participant-spectator matter 
and its development in terms of collision dynamics. Recently, participant- 
spectator matter has been found to be independent of the system size at 
energy of vanishing flow. Therefore, it acts as a barometer for the study of 
disappearance of flow and balance energy [TOl [TT] Section 2, describes the 
model briefly. In section 3, we present our results and section 4 summarizes 
the outcome. 

2 The model 

We simulate the nucleons within the framework of quantum molecular dy- 
namics (QMD) model. In the QMD model, [21 [3] each nucleon propagates 
under the influence of mutual interactions. The propagation is governed by 
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the classical equations of motion: 
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where H stands for the Hamiltonian which is given by: 
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The V i yrme , V^ uk , and Vf oul in Eq. (2) are, respectively, the Skyrme, 



Yukawa, and Coulomb potentials. 

3 Results and Discussion 

We here simulated the reactions of 40 Ca+ 40 Ca and 131 Xe+ 131 Xe, at different 
colliding geometries. This varies from very central (6 = fm) to peripheral 
one (b = Ri + R 2 ; Ri is the radius of either target or projectile). The incident 
energy is also varied between 200 and 400 MeV/nucleon. In the present study, 
we use a hard equation of state along with nucleon-nucleon cross section 
a = 40 mb. As discussed by many authors, the dynamics is insensitive 
towards the nuclear equation of state as well as towards nucleon-nucleon 
cross section. However, the above equation of state and cross section has 
been reported to reproduce the balance energy over wide range of colliding 
masses. [11] 

The present participant-spectator matter demonstration is based on the 
definition reported in our earlier recent publication. [TT] Here participant- 
spectator matter is defined in two different ways, (a) In the first definition, 
all nucleons experiencing at least one collision are labeled as participant mat- 
ter. The remaining nucleons are the part of spectator matter, (b) In the sec- 
ond definition, the above demarcation is based on the experimental method, 
where one uses different rapidity cuts to define participant and spectator 
matter. The rapidity of ith particle is defined as 



where E(i) and p z (i) are, respectively, the total energy and longitudinal 
momentum of iih particle. Now one can impose different cuts to study the 
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different participant-spectator matter. We shall here use the first definition 
to construct the participant-spectator matter. However, one should note that 
both these definitions have been reported to give same results. [TO] 

In Figs. [U and EJ we display the time evolution, respectively, in the 
spatial and momentum spaces for the central collision of 40 Ca+ 40 Ca at 200 
MeV/nucleon. The solid (open) circles denote the spectator (participant) 
matter. During initial stages of the reaction (till 10 fm/c), projectile and 
target are well separated both in the coordinate and momentum spaces. As 
a result, whole of the matter is in the form of spectator matter only. As the 
nuclei overlap, the density increases, that leads to more and more nucleon- 
nucleon binary collisions. This changes the matter from spectator to partic- 
ipant one. Due to the frequent nucleon-nucleon collisions, nucleons scatter 
in the transverse direction also. With the evolution of the reaction, partici- 
pant matter increases, leading to more nucleons scattering in the transverse 
direction. At the end of the reaction, we see only a few spectator nucleons. 

Looking at the momentum space (Fig. [2]), we see that the spectator nu- 
cleons have very little change in their initial velocity profile over the period 
of the reaction. On the other hand, participant matter leads the matter into 
midrapidity region. This drives the system into global equilibrium. Natu- 
rally, more the initial correlations are destroyed, more is the system close to 
the global equilibrium. Since, these figures are the snapshots of a single event, 
the picture could be different for different events. The matter is cold and 
fragmented at the final state. This may lead to the idea that compressed 
matter expands and due to Coulomb instabilities, this breaks into pieces. 
However, this argument has been questioned by one of us and collaborators 
who were able to detect the fragments at the time when matter is compressed 
and hot.[12j If one evolutes the reaction between time steps 10-30 fm/ c, one 
sees that the density rises from 0.9 to 1.7 between 10-20 fm/c, and falls by 
the same amount between 20-30 fm/ c. However, the average collision num- 
ber is not symmetric. The collision rate between 20-30 fm/c is higher than 
between 10-20 fm/ c. This points toward the compactness of the nucleons in 
space with slower velocities. 

In Figs. [3] and HJ we display the same evolution, but at colliding geome- 
tries of b/bmax = 0.2, 0.65, and 0.95, respectively. The solid (open) symbols 
represent spectator (participant) matter. The left panels of the figures 
are for the participant-spectator matter formed at the time of high density. 
The right panels are for the participant-spectator matter emerging at di- 
luted phase of the density. Note that at this state, the matter is cold and 



4 



40 



20 



- "Ca+^Ca 



- 



-20 - 



-40 



20 - 



S 
X 



-20 - 



-40 



20 - 



- 



-20 - 



-40 



0fm/ C ib/b ma =0.2 



E=200 MeV/nucleon 



+ 



(a) 



20 fm/c-L 



50 fm/c - 




1 1 — 

10 fm/c. 



<p ma *>/p„ =0.9 



<Ncoll> =0.01 



- 



# Spectator Matter 
(b) O Participant Matter' 
I I ' I I 



- -20 



-| 1 — 

30 fm/c. 



max , 

<P >/P„ 



=0.9 
=1.9 



<Ncoll> 
(d) . . avg . 

1 I I J I I I 

100 fm/c -| 

o 0)0 

oo o 

6> o o : o 
o qP&©6 o -I -20 



(0. 



o 



J i_ 



40 



20 



-40 

20 


- -20 
-40 



- 20 



- 



-40 



-40 



-20 



20 -40 -20 

z(fin) 



20 



40 



Figure 1: The time evolution of a single event of 40 Ca+ 40 Ca reaction at 200 
MeV/nucleon and b/b max = 0.2 in coordinate space. The solid (open) circles 
represent spectator (participant) matter, respectively. The < p max > / po 
and < Ncoll > avg are the maximum density and average number of nucleon- 
nucleon collisions averaged over large number of events. 



fragmented. The reaction time for the high densities is around 20 fm/c for 
b/bmax = 0.2, 0.65, and 0.95 as is also evident from Fig. [TJ Note that the time 
of maximal density varies linearly with the size of interacting nuclei. [5] [10] 
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Figure 2: Same as Fig. dj but in momentum space. 



From the figures, one notices a reduced intensity of the participant mat- 
ter with increase in the impact parameter. At central impact parameters, 
nucleons are well scattered in the momentum space. However, one can see 
two well separated momentum spaces at higher impact parameters. This 
indicates that momentum correlations remain preserved till the end of the 
reaction. Further, central collisions lead to significant portion as participant 
matter at high density. On the other hand, it is almost spectator matter at 
peripheral geometries. One may say that the reaction in central collisions is 
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Figure 3: The time evolution of a single event of 40 Ca+ 40 Ca reaction at 
200 MeV/nucleon. We use different impact parameters b/b max = 0.2 (upper 
panel), 0.65 (middle panel), and 0.95 (bottom panel). The solid (open) 
symbols, again, represent the spectator (participant) matter. The high, low, 
and med subscripts on the spectator/participant represent the intensity of the 
matter, being high, low or medium for a particular reaction. The displayed 
quantities like < p max > / p , < Ncoll > avg: and <R a > (defined in Eq. (4)) 
are averaged over hundreds of events. 
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Figure 4: Same as Fig. [31 but in momentum space. 



driven by the nucleon-nucleon binary collisions. The mutual two and three- 
body interactions dominate the physics at peripheral collision. Interestingly, 
though, the high density phase is nearly the same at b/b max = 0.2 and 0.65, 
the final stage evolutions is quite different. In the formal one, it is more of 
complete disassembly whereas it is a fire-ball picture in the later case. This 
demonstration also points toward a very important clue for multifragmenta- 
tion of these nuclei. Almost complete participant matter in central collisions 
does not allow any initial memory and correlations to survive. Whereas a 
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Figure 5: The trajectory of a single spectator nucleon and participant nu- 
cleon for the collision of 131 Xe+ 131 Xe at incident energy of 400 MeV/nucleon. 
The quantities in the boxes have superscripts and subscripts, representing, 
respectively, the number of collisions and time interval of these collisions. 
The shaded box is the interval for maximal density. 



complete spectator matter at the peripheral geometries preserves the initial 
correlations. Therefore, in both these situations, one should not expect in- 
termediate mass fragments (5 < A < 65) to be emitted. The formation of 
the intermediate mass fragments needs both the participant and spectator 
matter in a mild quantity. The semicentral collisions are, therefore, prefect 
for the study of emission of intermediate mass fragments. This has also been 
reported as rise and fall of intermediate mass fragments in several experimen- 
tal and theoretical studies. [131 E] Similar rise and fall also occurs with the 
change in the incident energies of the projectile. [12] Interestingly, baring the 
central collisions, no equilibrium is seen at semicental and peripheral geome- 
tries. Further from the Fig. HJ it is also evident that no substantial changes 
occur after the maximal high density in semicentral and peripheral collisions. 
However a nearly complete equilibrium can be seen in central collisions. 
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Figure 6: Same as Fig. but in momentum space. 



Let us now look for the changes in the phase-space with reference to 
the nucleon-nucleon collisions. In Figs. |5] and EJ we display the coordinate 
and momentum space, respectively, for a single spectator nucleon (denoted 
by the solid circle) and participant nucleon (denoted by the open circle). 
The participant nucleon has suffered a large number of binary collisions. 
The numeric values in the boxes represent the number of collisions faced 
by a nucleon. The subscripts to these values denote the time interval (in 
fm/ c) during which these collisions happened. From Fig. EJ one notices 
that the trajectory of spectator nucleon remains unchanged during the entire 
duration of the collision. The participant nucleon, however, shifts in the 
transverse direction after suffering frequent collisions. Similarly, from Fig. 
[6j momentum space of the spectator nucleon remains unchanged. Whereas 
in case of participant nucleon, ^-component of the velocity decreases and 
z-component of the velocity increases. The interesting case is the velocity 
profile between 19-20 fm/ c. Here, velocity is negative, leading to the reversed 
z-coordinate between this interval. This time also happens to be the time 
for the maximal density. We see that the participant matter generates also 
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Figure 7: The anisotropy ratio <R a > and relative momentum <Kr> for 
the case of 131 Xe+ 131 Xe at b/b max = 0.2. Here different lines represent the 
matter having different average collisions at different incident energies. The 
shaded area represents the time of high density > po for different reactions 
displayed here. 

the positive transverse momentum. 

In Fig. [71 we display the time evolution of anisotropy ratio <R a > and 
relative momentum <Kr> for different average collisions experienced in a 
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reaction. The <R a > ratio is defined as 



'vl + 



< Ra >= V ^ * ■ ( 4 ) 



This anisotropy ratio is an indicator of the global equilibrium of the system. 
This represents the equilibrium of the whole system and does not depend on 
the local positions. The full global equilibrium averaged over large number 
of events will correspond to <R a > = 1- The second quantity, namely the 
relative momentum <Kr> of two colliding Fermi spheres, is defined as 

< K R >=< \Pp(f,t) - P T (f,t)\/h >, (5) 

where 

E? 1 Pi(t)Pi(r,t) 

— * 

Here Pj and pj are the momentum and density of the jth particle and i 
stands for either projectile or target. The <Kr> is an indicator of the local 
equilibrium because it depends also on the local position r. As expected, 
the anisotropy ratio increases or remains almost constant whereas relative 
momentum decreases with increase in the number of collisions indicating 
better degree of global thermalization. The smaller value of (<K R >) also 
indicates toward the better thermalization of the matter. However, once 
nucleon-nucleon collisions exceeds a certain limit, the thermalization does 
not improve with further binary collisions. Rather the nucleons prefer trans- 
parency in the system. Naturally, the higher energy projectile has a lower 
value of (<R a >) at the start [conversely, the larger value of (<K R >)]. Fur- 
ther, as is also evident, (<R a >) ratio saturates as soon as high density phase 
is over. In other words, nucleon-nucleon collisions happening after high den- 
sity phase do not change the momentum space significantly. One also notices 
that in the absence of nucleon-nucleon collisions (in Vlasov mode), almost 
no change is seen in <R a >. 

In Fig. |HJ we show the normalized rapidity distribution, for a single 
spectator and participant nucleon. Although, the dN/dY will have the same 
value for both these particles, we have shifted the spectator and participant in 
y-axis for better clarity. Again, the values in the boxes represent the number 
of collisions suffered by a nucleon. The subscripts to these values represent 
the time interval (in fm/c). One notices that the participant nucleon, after 



12 



"i 1 r 



■ 131 Xe + 131 Xe 



i 1 r 



5 - 



^ 3 

7 



2 - 



1 - 



<p >=0.22 GeV/c 

T: 



"1 — 1 — I — 1 — I — 1 — 

E=400MeV/nucleon- 



^0 



b/b =0.2 

max 



••'•••■r^v- 



J i I i I i L 



• Spectator 
o Participant ' 



_i I i_ 



-2.0 -1.5 -1.0 -0.5 0.0 0.5 1.0 1.5 2.0 



Y / Y 

cm. beam 



Figure 8: The normalized rapidity distribution (dN/dY)„ orm as a function 
of Y c m /Y beam . Here a single spectator nucleon and participant nucleon at 
different time steps is taken. The shaded box is the interval for maximal 
density. 



suffering significant binary collisions, shifts from the projectile rapidity to 
midrapidity region. On the other hand, the spectator nucleon remains in the 
projectile rapidity region during the entire reaction. 
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4 Summary 



We have presented the simulations of heavy-ion collisions in terms of participant- 
spectator matter. We find that the participant-spectator matter depends 
crucially on the collision dynamics and on history of the nucleon. We see 
that the important changes in the momentum space are due to the binary 
collisions experienced during the phase of high density. This otherwise was 
not possible with the mean field alone. The collisions push the nucleons into 
midrapidity region that is responsible for the formation of participant mat- 
ter. This also indicates thermalization in heavy-ion collisions. Interestingly, 
large number of collisions, nevertheless, do not always guarantees a better 
degree of thermalization of the reaction. 
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